Topological characterization of periodically-driven quantum systems 
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Topological properties of physical systems can lead to robust behaviors that are insensitive to 
microscopic details. Such topologically robust phenomena are not limited to static systems but 
can also appear in driven quantum systems. In this paper, we show that the Floquet operators 
of periodically driven systems can be divided into topologically distinct (homotopy) classes, and 
give a simple physical interpretation of this classification in terms of the spectra of Floquet oper- 
ators. Using this picture, we provide an intuitive understanding of the well-known phenomenon of 
quantized adiabatic pumping. Systems whose Floquet operators belong to the trivial class simulate 
the dynamics generated by time-independent Hamiltonians, which can be topologically classified 
according to the schemes developed for static systems. We demonstrate these principles through 
an example of a periodically driven two-dimensional hexagonal lattice model which exhibits several 
topological phases. Remarkably, one of these phases supports chiral edge modes even though the 
bulk is topologically trivial. 



I. INTRODUCTION 

Following the discovery of the quantized Hall effect in 
198G 1 , there has been great excitement about the possi- 
bility of observing extremely robust, "topologically pro- 
tected" quantum phenomena in solid state systems^. 
This excitement has been redoubled in recent years with 
the discovery of new classes of materials called "topolog- 
ical insulators" in two and three dimensional systems^—. 
A common feature linking all of these systems, which 
from a traditional point of view appear to be mun- 
dane band insulators, is the fact that their ground state 
wave functions feature internal structures characterized 
by non-zero values of integer topological invariants, which 
distinguish them from conventional, trivial systems. The 
existence of such internal structures leads to the appear- 
ance of robust gapless edge modes wherever a non-trivial 
material has an interface with a trivial one, such as the 
vacuum. Due to their topological origin, these modes can 
not be localized or destroyed by a wide range of perturba- 
tions, and give rise to many interesting, robust phenom- 
ena. Recently, several groups have proposed a compre- 
hensive scheme for classifying all possible types of such 
topological phases which can arise in band insulating and 
superconducting systems^— . 

Meanwhile, a variety of robust topologically-protected 
phenomena have also been found to occur in the dy- 
namics of driven quantum systems with time-dependent 
Hamiltonians. Such phenomena can roughly be divided 
into two broad classes. First, there is a class of phenom- 
ena displaying quantized adiabatic pumpingii~— , which 
can be conceptually traced to Thouless' original proposal 
of quantized adiabatic transpor t 14 i 15 . Second, many 
groups have studied the possibility of effectively simulat- 
ing the behavior of the topologically non-trivial (static) 
materials described above by applying periodic driving 
fields to artificial^— or condensed matter systems^—. 
In this paper we provide a unifying scheme for character- 
izing these various types of topological phenomena which 
occur in the dynamics of periodically driven systems. 



Furthermore, from this unified view, we find that these 
two classes of phenomena are not unrelated, and can in 
fact both be realized within a single periodically-driven 
graphene-like system. 

Unlike the case of static (non-driven) systems, 
periodically-driven systems do not have well-defined 
ground states which can be used for classification. In- 
stead, we classify driven systems in terms of the topo- 
logical properties of their corresponding "Floquet opera- 
tors," i.e. their time evolution operators acting over one 
full period of the drive, T. Each eigenstate of the Floquet 
operator, called a Floquet state, accumulates a phase 
<p over one period of the driving. Accordingly, to each 
Floquet state, we associate a "quasi-energy" e — <p/T, 
which is the average phase accumulated per unit time. 
For many purposes, the Floquet states and their associ- 
ated quasi-energies can be regarded in an analogous way 
to the eigenstates and corresponding energies of a static 
system. However, because the quasi-energy is defined as 
a phase variable, it is periodic with period 2tt/T. This pe- 
riodicity introduces an additional topological structure, 
associated with the winding of quasi-energy, which has 
no analogue in static systems. As we will discuss below, 
this property allows a number of interesting phenomena 
to occur in driven systems, such as quantized pumping, 
and even the existence of chiral edge modes in a two- 
dimensional system with topologically-trivial bulk bands 
(i.e. with all Chern numbers equal to zero). 

The mathematical structure which captures these 
quasi-energy winding-related phenomena is that of 
the homotopy groups of unitary (Floquet) operators. 
Throughout the main text, we will focus primarily on 
systems with discrete (lattice) translational symmetry in 
d-dimensions, which are subjected to spatially homoge- 
nous, periodic, time-dependent driving. In this case, we 
obtain simple expressions for topological invariants asso- 
ciated with the first and third homotopy groups, writ- 
ten in terms of Floquet operators parameterized by the 
conserved crystal momentum k. The case of disordered 
systems is discussed in terms of a "twisted boundary con- 
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dition" approach-^ in Appendix [S] 

When the homotopy group classification returns a triv- 
ial result, i.e. for systems without quasi-energy winding, 
the Floquet operator can be expressed in terms of a lo- 
cal effective Hamiltonian H e s through U(T) = e ~ lHoltT . 
Here one can view the dynamics as a stroboscopic simu- 
lation of the dynamics of a static system governed by the 
Hamiltonian H c h . In this case, the topological character- 
istics of the driven system can be classified in terms of the 
symmetry and dimensionality according to the schemes 
laid out for static topological insulators and supercon- 
ductors in RefsisMTol. Just as in the static case, systems 
characterized by non-zero values of the topological invari- 
ants in these schemes possess robust chiral edge modes 
at interfaces with "trivial" regions^. 

The characterization of periodically driven systems in 
terms of the topological structure of Floquet operators 
constitutes the major result of this paper. This approach 
provides a natural description of topologically-quantized 
pumping, and reveals a simple and intuitive picture in 
which to understand this phenomenon. Furthermore, the 
general mathematical structure provides a guide for iden- 
tifying new classes of topologically-protected behavior 
which may arise in driven systems. The intuition gained 
from this approach may help open new avenues in which 
to explore topologically-robust behaviors in artificial sys- 
tems as well as in natural materials, where experimental 
realizations are already feasible with current technology. 

The paper is structured as follows. In Section UH we 
define the class of driven systems that we will consider, 
and review the basic language of Floquet theory. Then 
in Section Mil we present the classification of driven sys- 
tems in terms of the homotopy groups of their associated 
Floquet operators. Here we show that the topological 
invariant V\ associated with the first homotopy group of 
the Floquet operator directly measures the winding of 
the quasi-energy spectrum, and describe how it leads to 
the quantized average displacement. In Section lTVl we de- 
scribe the second scheme of topological classification of 
driven systems based on "effective Hamiltonians," which 
applies in the absence of quasi-energy winding. Here we 
provide a detailed discussion of the meaning of effective 
time reversal and particle-hole symmetries for driven sys- 
tems. In Section [V] we demonstrate the appearance of 
both topological structures described in Section [Hi] and 
IIVI in the dynamics within a two dimensional tight bind- 
ing model on a hexagonal lattice with time dependent 
hopping amplitudes. Strikingly, over a range of param- 
eters, this model exhibits a phase which features chiral 
edge modes in the absence of topologically non-trivial 
bulk quasi-energy bands. We conclude in Section |VT1 



II. FLOQUET THEORY FRAMEWORK 

In this section, we describe the Floquet theory frame- 
work that will be used to analyze the dynamics of driven 
systems in this paper. We focus on the driven dy- 
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FIG. 1: Typical schematic instantaneous band structure con- 
sidered in this work, illustrated for a one-dimensional system 
with crystal momentum k. We consider the case where m 
low-lying bands are separated from higher bands by a band 
gap Eg at all times t during the periodic driving cycle. We 
assume that the frequency of driving is much smaller than the 
band gap Eg at each instant, such that excitations out of the 
m-fold low-lying band subspace can be neglected. 



namics of a quantum system of non-interacting particles 
in a d-dimensional lattice, subjected to a periodically- 
varying, time-dependent Hamiltonian H(t), which satis- 
fies H(t + T) = H{t). Here T is the period of the driving 
cycle. The behavior is analyzed in terms of the evolution 
operator of the system over one full period of the driving, 
U(T), defined as 



U{T) = Te- i f° H( - t)dt , 



(1) 



where T is the time-ordering operator. In the theory of 
periodically-driven systems^, U(T) is called the Floquet 
operator. Note that while U(T) implicitly depends on the 
starting point of the interval, here taken to be t = 0, the 
topological classification described below is independent 
of the choice of the starting point. 

Throughout the main text, we assume that the sys- 
tem possesses discrete (lattice) translational symmetry 
at each instant of time. However, topological classifi- 
cation can still be applied even when the system lacks 
translational invariance. The extension to this more gen- 
eral case is discussed in Appendix [S] In the presence 
of translational invariance, the crystal momentum k is a 
conserved quantity of the dynamics. Therefore, at each 
time t, the Hamiltonian H(t) can be characterized by a 
band structure as shown schematically in FigHJ 

Suppose that at time t = 0, the state of the system 
lies entirely within the subspace spanned by the m low- 
est bands of the instantaneous Hamiltonian H(0). In 
this work we consider the situation where the driving 
always returns the system to this subspace after each 
full period, without causing excitations to other bands. 
This condition is guaranteed, for example, when higher 
energy bands of H(t) are well separated from lower en- 
ergy bands by a large energy gap E g (see FigH]) and the 
evolution is adiabatic with respect to E g at all times. 
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In this situation, the evolution operator for one full pe- 
riod of driving is described by the set of to x to matrices 
{[/ k (T)}, labelled by the value of the crystal momentum 
k. In terms of these operators, the Floquet operator, 
Eq-CQ), projected onto the low-energy subspace, is given 
by U(T) = Ek C/ k( T ) ® Pk, where P k is the projector 
onto the TO-dimensional subspace of low-energy states 
with crystal momentum k. Note that for intermediate 
times < t < T, even in the fully adiabatic regime, the 
evolution operator U^{t) will in general include mixing of 
higher bands of the initial Hamiltonian H(0) and there- 
fore cannot be represented in the basis of the m lowest 
bands of H(0). The restriction to to dimensions for H(0) 
applies only for complete cycles of the periodic driving. 

In analogy with the energy associated with each eigen- 
state of the Hamiltonian of a time-independent system, 
we associate a quasi-energy e with each eigenstate \<p) of 
the Floquet operator of the driven system. The quasi- 
energy is defined through the phase accumulated by the 
"Floquet state" \tp) over each full period of driving as 



a) Quasi-energy (l/T) b) Quasi-energy (1/T) 



U{T)W) 
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\tp). Note that, unlike real energies 



which are uniquely determined, quasi-energies are only 
uniquely defined up to integer multiples of 2-k/T. This 
property will be crucial to the analysis below. For sys- 
tems with discrete translational symmetry, we index the 
quasi-energies by the crystal momentum k, and define a 
"quasi-energy band structure" {£k,a}, where a is a band 
index. In the next section, we will see that the "winding" 
property of quasi-energy allows for quasi-energy band 
structures, and hence behaviors, which are qualitatively 
different from those found in static systems. 



III. HOMOTOPY GROUP CLASSIFICATION 
OF EVOLUTION OPERATORS 

In this section, we formulate the topological classifi- 
cation of periodically driven systems in terms of the ho- 
motopy groups of the evolution operators {J7k(T)}. In 
subsection IIII A[ we focus on the topological invariant 
v-i associated with the first homotopy group of evolution 
operators. We illustrate many general properties of the 
topology in evolution operators through the study of an 
example in one dimension. Then in subsection IIII Bl we 
discuss generalizations to higher dimensions. 



A. Topological invariant v\ and quasi-energy 
winding 

In one spatial dimension, the set of operators {Uk(T)} 
defines a map from the Brillouin zone — 7r < k < n to the 
space of m x m unitary matrices. Due to the periodicity of 
crystal momentum, the Brillouin zone is a circle, and this 
map traces out a closed loop in the space ofmxm uni- 
tary matrices. This loop can be characterized in terms of 
a homotopy class, which identifies all such loops that can 
be smoothly deformed into one another. These homo- 
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FIG. 2: a) Quasi-energy spectrum of a spin-1/2 particle in 
a time- and spin-dependent potential, displaying non-trivial 
topology characterized by v\ = 1 (see text). The quasi-energy 
is plotted in units of 1/T. The topological invariant v\ counts 
the total winding of the bands in the quasi-energy direction 
as crystal momentum k is taken around the Brillouin zone, 
b) The total winding number v\ is unchanged when coupling 
between up and down spin states is introduced. 



topy classes are indexed by an integer-valued topological 
invariant v\ that captures the topology (or "winding" ) of 
the corresponding maps. 

In the presence of translational symmetry, the invari- 
ant v\ is defined as: 



1 

2^ 



dkTr 



U k {T)-Hd k U k {T) 



(2) 



where k is integrated over the first Brillouin zone, and 
the trace is taken over the TO-dimensional internal space 
of the Bloch wave functions^. We emphasize that a 
map characterized by the topological invariant v\ can 
be defined in a system with spatial dimension d > 1 by 
choosing k in Eq.Q to parametrize some closed, non- 
contractible loop in the d-dimensional Brillouin zone. 
For example, consider the case of a two-dimensional sys- 
tem on a square lattice with crystal-momentum compo- 
nents k x and k y . Here the Brillouin zone is a torus, 
and we can construct two invariants v\ x and v\ y as- 
sociated with two distinct loops on the torus, k x = 
constant, — 7r < k y < n, and — ir < k x < n, k y = constant: 

v Xx = ± J^dkyT^UkiTyHd^UkiT)} for a fixed k x , 
and similarly for v\ y . Note that, due to continuity, v\ x 
(y\ y ) is independent of the value of k x (ky). Hence this 
first homotopy group classification is not restricted to 
systems in one spatial dimension. 

The meaning of the topological invariant v\ is best 
illustrated through a simple example. Consider a spin- 
1/2 particle in a deep one-dimensional, spin-dependent 
lattice potential. Here and throughout the paper, we take 
the lattice constant to be 1. The potential for the spin- 
down state is time- independent, while the potential for 
the spin-up state moves slowly to the right, such that over 
the time T it shifts by exactly s times the lattice constant, 
where s is an integer. If we assume that the movement 
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is adiabatic, then the Floquet operator restricted to the 
space of the lowest band of initial Hamiltonian H(0) is 
given by U(T) = J2 X \ X + s )( x \ ® p t + 1 ® P U where 
■Per = I^X "! is a projector onto the spin state a =t,l, 
and \x) is a state localized in a single well of the (deep) 
lattice potential, with x an integer labeling the unit cell. 
In the crystal momentum basis, we have Uk(T) — e~ lsk ® 
Pf + 1 O Pj,. This evolution is indeed characterized by 
the non-trivial value fi = s, as can be checked explicitly 
using Eq.©. 

We can obtain an intuitive understanding of v\ by eval- 
uating Eq.@ in the basis of Floquet states. This analysis 
gives 

where {£k,a} are the quasi-energy bands labeled by the 
band index a 29 . The integral in Eq. <j3j> counts the total 
winding number of the quasi-energy bands as k runs over 
the first Brillouin zone. This picture of v\ is particularly 
appealing since it allows the non-trivial topology to be 
identified through simple inspection of the quasi-energy 
spectrum. The quasi-energy spectrum for the case s = 1 
is plotted in Figj2£i. Note that such winding is only pos- 
sible due to the periodicity of quasi-energy for a driven 
system; therefore v\ = for any local, static (non-driven) 
system. 

An important physical manifestation of the topology 
captured by v\ is revealed by Eq.((3]). The factor dsk, a /dk 
plays the role of a group velocity for the Floquet band a. 
Due to the periodicity of Ek, a in k and the periodic nature 
of quasi-energy, the average group velocity is quantized: 
the average slope of e^.a must take a value which yields 
an integer number of windings over the Brillouin zone. 
Averaged over all k, the displacement Ax = (dsk, a /dk) ■ 
T after a full period T is therefore quantized. In the 
simple example above, in which the particle is initially 
in a uniform superposition of states with all values of k, 
this quantization means that over one period the average 
position of a particle in the spin-up band is shifted by 
exactly s unit cells to the right, while the average position 
of a particle in the spin-down band does not change. 

For the case of adiabatic evolution of a filled-band 
fermionic system, it can be shown that v\ is directly re- 
lated to the charge current integrated over one period. 
Therefore, in this case, the quantization of v\ implies 
the quantization of pumped charge first identified by 
Thouless 14 . In Appendix IB1 we explicitly show this rela- 
tion. Moreover, as is implied by the relation between v\ 
and quantized charge pumping, v\ can be directly related 
to the first Chern number of the Bloch wavefunctions of 
the filled bands. We demonstrate the general relation 
between the topology captured by homotopy groups of 
Floquet operators and Chern numbers in Appendix [C] 

While the homotopy-class-based topological character- 
ization of periodically driven systems given in this pa- 
per is closely related to the Chern-class-based charac- 
terization of adiabatic pumping developed in previous 



work a 10 ! 1 , we emphasize that our approach provides an 
intuitive framework which naturally leads the way to gen- 
eralizations to systems in other dimensions or with ad- 
ditional symmetries (see Section fill Bp . Moreover, as il- 
lustrated by an example in Section [Vj the picture based 
on the winding of quasi-energy makes it possible to iden- 
tify non-trivial topological behavior of driven systems in 
a simple, direct manner. 

The integer-quantization of v\ implies that its value is 
robust against various perturbations. For instance, the 
value of v\ is insensitive to continuous deformations of 
the quasi-energy spectrum, and to mixing with other, 
topologically trivial, quasi-energy bands. For the model 
described above with quasi-energy band structure shown 
in Figj^, such mixing can be introduced, for example, 
by applying a spin rotation pulse during the dynamics. 
The quasi-energy spectrum in the presence of such mix- 
ing is shown schematically in FigJSJa. Note that, even in 
this case, Eq.Q for the net winding number still gives 
i/i = l; hybridization of Floquet bands with trivial and 
non-trivial topology does not change the total winding 
number. 

In the presence of disorder, translational invariance is 
destroyed. The crystal momentum k is not a good quan- 
tum number, and the definition of v\ in Eq.@ cannot be 
used. Nonetheless, the topologically protected phenom- 
ena are expected to be robust against weak perturba- 
tions, including disorder. Indeed, in Appendix IA1 we use 
the method of twisted boundary conditions^ to define 
a generalization of the topological invariant V\ , which is 
valid even in the presence of disorder. This generalized 
invariant can also be used to classify the dynamics of peri- 
odically driven, gapped, interacting many-body systems, 
such as filled-band states of fermions or Mott insulating 
states of bosons^ 2 -. 



B. Higher dimensional systems 

In the previous subsection, we saw that the phe- 
nomenon of quantized adiabatic pumping in one dimen- 
sion acquires a simple and intuitive explanation when an- 
alyzed in terms of the topological structure of Floquet op- 
erators. We now discuss how this elegant mathematical 
framework can guide us in searching for generalizations of 
topological pumping in higher dimensional systems. The 
key observation is that the Floquet operators of period- 
ically driven systems with dimension greater than one 
can be characterized by topologies associated with ho- 
motopy groups beyond the fundamental group. Here we 
describe the general situation in higher dimensions, and 
give an explicit expression for the topological invariant 
associated with the third homotopy group, which can be 
relevant in three-dimensional systems. 

For a system in d dimensions, the first Brillouin zone 
is topologically equivalent to a d-dimensional torus. By 
allowing k to parametrize an ^-dimensional toroidal sec- 
tion of the Brillouin zone, with £ < d, we can use the 
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set of (Floquet) evolution operators {U^T)} to define a 
map from the ^-dimensional torus to the space ofmxm 
unitary matrices, U(m). Such maps can be classified us- 
ing the structure of the £-th homotopy group of U(m), 
denoted by 7i^[U(m)]. Noting that the homotopy groups 
are given by 7r<2[U(m)] = Z for odd d and 7Td[U(m)] = f 
for even d, with m > (d + l)/2— , we see that a new 
type of topological invariant appears in each odd dimen- 
sion. While the d-th homotopy group of generic unitary 
matrices is trivial when d is even, the Floquet opera- 
tors of even-dimensional systems can still possess non- 
trivial topological structures either under lower homo- 
topy groups (e.g. 7Ti[U(m)] for a system in d — 2), or if 
the space of allowed evolution operators is restricted by 
additional symmetries of the system. 

For systems in three dimensions, d = 3, we consider 
the topological invariant vj, associated with the third ho- 
motopy group: 



Tr 



d A k 
24tt 2 

-l 



(4) 



U?d k JJv) {U^d kfi U k ) (U^d^U*) 



where £ Q ^ 7 is the anti-symmetric tensor, a, (3, 7 = x, y, z, 
and k is integrated over the first Brillouin zone. Note 
that vz can only be non-vanishing if the dimension of 
Uk(T) is larger than 1 so that multiplication can be non- 
commutative. The search for physical realizations and 
the manifestations of this topological invariant will be 
interesting subjects for future work. 

The phenomena associated with non-trivial topology of 
Floquet operators under the homotopy groups are unique 
to periodically-driven systems. For a system governed 
by a static Hamiltonian Hs, the eigenstates of the Flo- 
quet operator Us(T) — e ~ zHsT and the associated quasi- 
energies for any choice of T coincide with those of Hs- 
Because the energy eigenvalues E^ n of any local Hamil- 
tonian Hs are smooth, periodic functions of k, the corre- 
sponding quasi-energy eigenvalues £k,n cannot wind and 
there must exist a gap in the quasi-energy spectrum. It 
can be shown that a "gapless" quasi-energy spectrum, 
in which there is at least one eigenstate for every quasi- 
energy — ir/T < e < x/T, is a necessary requirement for 
non-trivial topology under any homotopy group^. An 
example of such a gapless spectrum is shown in Fig. [5] 
Therefore the Floquet operator of any local, static system 
is trivial under all homotopy groups. For a periodically- 
driven system with driving period T, however, the quasi- 
energy spectrum associated with the evolution operator 
U (T) need not have a gap. Hence, as demonstrated by 
the example above, periodic driving can produce evolu- 
tion operators with non-trivial topology under the homo- 
topy groups. 



IV. TOPOLOGICAL INVARIANTS OF GAPPED 
EFFECTIVE HAMILTONIANS 

So far, we classified the dynamics of periodically driven 
systems with discrete translational symmetry through 
the generalized "windings" of their Floquet evolution op- 
erators as identified by the homotopy groups of unitary 
matrices. However, even when there is a gap in a sys- 
tem's quasi-energy spectrum and its Floquet operator is 
thus trivial under all homotopy groups, its dynamics may 
possess other topological characteristics associated with 
an effective Hamiltonian H c h: 



U(T) = Te~ l $« H{t)dt ee e~ iHettT . 



(5) 



When viewed in this way, evolution under the period- 
ically varying Hamiltonian H(t) stroboscopically simu- 
lates the evolution of a system with a static Hamiltonian 
H c s at integer multiples of the driving period T . 

When H e g has a band structure with a band gap, the 
previously introduced classification schemes^— for static 
insulators and superconductors can be directly applied 
to H e ff^. In these schemes, systems are categorized into 
distinct topological classes based on dimensionality and 
on the presence or absence of the (discrete) time-reversal 
and particle-hole symmetries. In a periodically-driven 
system, these symmetries must be considered at the level 
of the effective Hamiltonian, or equivalently, of the Flo- 
quet operator U(T), as we now describe. 

At the level of the Floquet operator U(T), time- 
reversal symmetry is defined by the existence of a unitary 
operator Q which has the action QH* S = H c g, or 



QU{T)*Q) = [/(T)t. 



(6) 



Here H* s is the complex conjugate of H c g . The operator 
Q satisfies QQ* = ±1, where the sign, + or — , distin- 
guishes the cases with integer or odd-half-intcger spin, 
respectively. A sufficient, although not necessary, condi- 
tion for the presence of time-reversal symmetry is satis- 
fied if the time-dependent Hamiltonian possess a special 
point to such that 



QH(t + t )*Qi = H(-t + t ). 



(7) 



We give a detailed proof of this statement in Appendix 
ID] Note that condition in Eq. ([7]) can be satisfied even if 
the instantaneous Hamiltonian H(t) is not time-reversal 
invariant, such as when the system is subjected to mag- 
netic fields. Conversely, even if H(t) instantaneously has 
time-reversal symmetry for each t, the condition in Eq. 
(|7|) is not necessarily satisfied, and H e g need not be time 
reversal invariant. This latter situation arises in the ex- 
ample which will be explored in SecfVl 

Similarly, particle-hole symmetry is defined at the level 
of the Floquet operator by the existence of a local unitary 
operator V satisfying 



VU(T)*V' f = U(T) 



(8) 
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with VP* = ±1. Any Bogoliubov-de Gennes Hamilto- 
nian describing the dynamics of superconducting quasi- 
particles possesses this symmetry (with W* = 1), even 
if a time-dependent perturbation is added. 

Many results from the study of topological properties 
in static systems can be directly translated to period- 
ically driven systems, with the understanding that the 
topologically-protected phenomena apply to the behav- 
iors of Floquet states. For example, the edge of a static 
two dimensional system characterized by a non-zero 
Chern number is known to host chiral edge states. In Sec- 
tion |V1 we study dynamics in a driven two-dimensional 
tight-binding system, and show that driving can induce 
non-zero Chern numbers in the bands of the effective 
Hamiltonian. The edge of such a system hosts chiral Flo- 
quet edge states which propagate unidirectionally when 
viewed at integer periods of the driving. Analogously, the 
edge of a system with a time-reversal-invariant effective 
Hamiltonian can host helical Floquet edge states^. 




FIG. 3: Hexagonal lattice structure. Sublattice A is marked 
with filled circles and sublattice B is marked with open circles. 
Ji for i = 1, 2, 3 represent the hopping amplitudes between the 
sites. 



nT< t <nT+T/3 



V. DYNAMICALLY-INDUCED TOPOLOGICAL 
PHASES IN A HEXAGONAL LATTICE 

In this section, we study single-particle dynamics in a 
two dimensional hexagonal lattice tight-binding model, 
where the hopping amplitudes are varied in a spatially 
uniform, but time-dependent, cyclic fashion. Here we 
choose the hexagonal lattice but expect that similar 
physics can be realized in other lattices as well. In differ- 
ent parameter regimes, the system can support topologi- 
cal phases of either the Floquet operator homotopy type, 
or of the effective Hamiltonian type. For weak driving, 
there is no winding of the quasi-energy bands {v\ = 0), 
but the driving induces non-zero first Chern numbers in 
the two bands of the effective Hamiltonian. For larger 
driving amplitudes, the Chern numbers associated with 
each of the two bands become zero. While the bulk topol- 
ogy given by the Chern number is trivial, the invariant 
V\ becomes non-zero in a system with edges. 

We start from a tight-binding model Hamiltonian on a 
hexagonal lattice: 



cub 



H = E ( C U <b ) B(k) 

k 

H(k) = — Ji(t) (cos(bi • k)a x + sin(b,; • k)<7 y ) 



(9) 



where A and B label the two sublattices of the hexag- 
onal lattice, {Ji} are the hopping amplitudes from B 
sites to the neighboring A sites in the three directions 
i = 1,2,3 (see Fig. [3]), and {b^} are the vectors given 
by bi = (-1/2.V3/2), b 2 = (-1/2,-^3/2) and b 3 = 
(1,0). Here is the creation operator for a Bloch state 
with crystal momentum k on sublattice a = A 7 B. Ex- 
pressed in terms of local (site-specific) creation operators 

{4 4 ,J: we have c l, a = T^Ex, 4,, Q e lk ' x % where N is 
the number of unit cells in the system. Here and in the 





nT+T/3< t <nT+2T/3 



m 

nT+2T/3<t<(n+l)T 



FIG. 4: Driving cycle considered in the text, in which the 
hopping amplitudes {Ji} are varied in a cyclic fashion. Here 
we consider only A > 1, where the hopping amplitude along 
one of the three bond types is uniformly increased relative to 
the other two during each stage of the cycle. 



following, the distance between any site and its nearest 
neighboring sites in the hexagonal lattice is taken to be 
1. 

We consider a driving protocol where the hopping am- 
plitudes {Ji(t)} are modulated cyclically in time accord- 
ing to (see Figj4j): 

1. Ji = AJ; J 2 , J 3 = J for nT < t < nT + T/3 

2. J 2 = AJ; Ji, J 3 = J for nT + T/3 < t < nT + 2T/3 

3. J 3 = AJ; J x , J 2 = J for nT + 2T/3 <t<nT + T. 
In the following, we consider the cases A > 1. Notice 



7 



that when A = 1, the Hamiltonian describes hopping on 
a time-independent hexagonal lattice with uniform hop- 
ping amplitudes. Just as in graphene, in this case the 
spectrum is gapless, with two inequivalent Dirac points. 
We focus on the high frequency limit, J/cj <C 1, where 
uj = 2ir/T is the angular frequency of the drive. In this 
limit, all states within the two bands are far from res- 
onance, and the graphene-like band structure provides 
a good zeroth-order approximation to the quasi-energy 
spectrum. 

We are interested in the topological properties of the 
dynamics when A is increased from 1. For A > 1, we 
obtain the effective Hamiltonian through 

e -iTHctf = e -iH 3 T/S e -iH 2 T/3 e -iHiT/3 ^ 

where Hj is the Hamiltonian during the time nT + (J — 
l)T/3 < t < nT + jT/3 (see protocol above). For A > 1, 
gaps open at the Dirac points in the quasi-energy spec- 
trum of H c ff. For small (A — 1) <C 1, the size of the gap 
is given by -\/3[(A — 1) J] 2 T. Because H e g is a two-band 
Hamiltonian, we can write^ 

ff eff (k) = e(k) n(k) • ct, (11) 

where er = (a x , <r y , <r 2 ) is a vector of Pauli matrices acting 
in the sublattice space, and n(k) is a three dimensional 
unit vector. When e(k) ^ 0, n/T for all k, each quasi- 
energy band of the effective Hamiltonian can be charac- 
terized topologically by the first Chern number. The first 
Chern numbers for the two bands are given by^: 

C± = ^ / n {8 k n x d ky n)d 2 k. (12) 

Here, the k integration is taken over the first Brillouin 
zone. 

We numerically calculate the Chern numbers C± from 
Eq. (fT2"j) using a fixed ratio J/oj — J| between the hopping 
strength J and the frequency lu. We find C± = ±1 for 
1 < A < A c , where A c ~ 3.3 is the critical coupling where 
a topological phase transition occurs. At A = A c , the 
quasi-energy gap at e = ±ir/T closes for the states at 
k = 0. For A > A c , the Chern numbers of both bands 
become zero until another phase transition point A = 
X' c , where X' c w 8.7 is reached^, 4 . Thus we see that the 
time-dependent driving within a hexagonal lattice tight- 
binding model can induce non-zero Chern numbers in the 
effective Hamiltonian. 

A consequence of non-zero Chern number in a sys- 
tem with boundaries is the existence of chiral Floquet 
edge states. For example, a system with a "strip" ge- 
ometry (see Fig. [5^) is expected to have two counter- 
propagating chiral edge modes localized on the upper and 
lower boundaries. In the following, we consider the strip 
geometry with edges of armchair type. 

In Fig. |H1 we plot the quasi-energy spectrum of the 
driven strip as a function of the crystal momentum par- 
allel to the strip, for several values of A. When A = 1, 




-ti/3 tt/3 



FIG. 5: a) The strip geometry considered in the text, with 
finite extent in the y direction, and armchair edges running 
parallel to i. b) Spectrum of the (static) armchair ribbon 
shown in panel a) with a width of 100 lattice sites in the y 
direction. Energy is plotted in units of J. 



the Hamiltonian is independent of time, and the quasi- 
energy spectrum is identical to the true energy spectrum 
of a graphene ribbon with armchair edges. For A = 3, 
chiral edge states are visible within the gap centered at 
quasi-energy e = 0. The states are color-coded such that 
states which are localized at the upper edge are shown in 
red (dark color) and those at the lower edge are shown 
in green (light color). Although the results are shown for 
an armchair- type edge, analogous edge states are formed 
for other edge types^ 

At A = A c , the quasi-energy gap closes at e — ±ir/T. 
Since the gap centered at e — remains open at this 
point, the chiral edge states cannot be removed as A is 
increased through A c . Remarkably, the system supports 
topologically-protected chiral edge states for A > A c , de- 
spite the fact that the Chern numbers of both bulk bands 
are zero. The chiral edge states in the quasi-energy spec- 
trum for A = 4 > A c are shown in FigJBJi). Isolating 
the two edges, note that each edge state has a "winding" 
structure which resembles that shown in Fig. [2^,). This 
suggests that these edge modes may be characterized in 
terms of the topological invariant v\ studied in Sec. IHI Al 

In order to understand the origin of these features, it 
is helpful to consider the limit J — > 0, while keeping 
AJT/3 = 7r/2 non-zero. In this limit, during each of the 
three stages of the cycle, only sites connected by the bold 
lines in Fig. U are coupled. The amplitude of hopping 
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FIG. 6: Quasi-energy spectrum of -ff e ff in a finite strip of 
width 100 lattice sites in the y direction with armchair edges 
along the x direction. Here we choose JT = n/16, and 
A = 1,3,3.3,4. For A = 3, the two bands of H e f; are char- 
acterized by non-zero Chern numbers, which are manifested 
in the presence of chiral edge modes. The lines correspond- 
ing to these modes are colored such that the mode localized 
on the upper (lower) edge is shown in red (green). Close to 
A = 3.3, the gap of H e g at quasi-energy n/T closes, and then 
opens again as A is increased further. At A = 4, the Chern 
numbers associated with each of the bands are zero, yet the 
nanoribbon clearly still supports chiral edge states. 



between the two lattice sites of a dimer is chosen such 
that after a time period of T/3, a particle is transferred 
with certainty from one lattice site to its neighbor. The 
dynamics in this limit is depicted in Fig. [7J where the 
red and green solid lines show the propagation of parti- 
cles along the upper and lower edges, respectively, and 
the blue dotted line shows the propagation of particles in 
the bulk. Note that all the Floquet states in the bulk arc 
localized: a particle starting at any site in the bulk comes 
back to the starting site after two complete driving pe- 
riods. On the other hand, a particle which starts at the 
upper (lower) edge on sublattice B {A) propagates unidi- 
rectionally along the edge to the left (right). Therefore, 
the bulk and edge states are clearly separated. 

In the limit described above, it is straightforward to 
confirm that the Floquet evolution operators {Uk„.(T)} 
restricted to the space spanned by sites on sublattice B 
{A) of the upper (lower) edge, yield v\ = 1(— 1). These 
chiral Floquet edge states cannot be removed unless the 
bulk quasi-energy gap closes. Thus, for any other values 
of J and A which can be reached from this limit without 
closing the bulk gap, the existence of chiral Floquet edge 
states is guaranteed. This is indeed the case for A = 4 
and JT — tt/16, whose quasi-energy spectrum is shown 
in FiglHl 



2T 




FIG. 7: Dynamics under the driving cycle shown in Fig[4] 
in the limit J ->■ 0, A -> oo with AJ(T/3) = tt/2. The red 
(green) solid line shows the unidirectional propagation of a 
particle initially localized on the B (A) sublattice along the 
upper (lower) edge. As indicated by the blue dotted line, all 
bulk states are localized. A particle initially localized at any 
site in the bulk returns to its original position after a time 
2T. 



Recently, there have been several proposals to induce 
topological phases with AC electromagnetic fields^— . 
In spirit, these models are similar to ours in the small 
(A — 1) limit, where the bulk bands possess non-trivial 
topology (in our case and in Re fl24l ch aracterized by the 
first Chern number, and in Rcfs. 25.26 characterized by a 
Zi invariant). In particular, all of these models support 
robust chiral or helical Floquet edge states. However, the 
existence of chiral edge modes without non-trivial topo- 
logical bulk quasi-energy bands is a unique feature of 
our model, which displays a non-trivial winding of quasi- 
energy. 



VI. DISCUSSION 

We have presented a unified framework for the topo- 
logical classification of periodically-driven systems. Such 
systems are characterized in terms of their Floquet oper- 
ators, i.e. their unitary evolution operators acting over a 
full period of driving. Non-trivial topology can arise on 
two levels. First, we characterize the topological struc- 
ture of the Floquet operator in terms of the homotopy 
groups of unitary matrices. Here, topologically distinct 
classes of dynamics are identified by the values of (gen- 
eralized) "winding numbers," v\ and 1/3 for the first and 
third homotopy groups. Explicit expressions for v\ and 
V3 are given for systems with discrete translational sym- 
metry in Eqs.([2|) and ([5]), and in terms of a more gen- 
eral twisted boundary condition formulation in Eqs. (|A2p 
and (IA4|) . Within this framework, the well-known phe- 
nomenon of quantized charge transport^ obtains a nat- 
ural and intuitive explanation in terms of the winding 
of quasi-energy. Systems with topologically non-trivial 



9 



Floquet operators in this homotopy sense can exhibit 
interesting features, such as a chiral dispersion in one- 
dimension, or protected edge states in two-dimensional 
systems with topologically trivial bulk band structures. 
Such phenomena cannot be realized in systems governed 
by local, static (time-independent) Hamiltonians. 

When the Floquet operator is "trivial" in this ho- 
motopy group sense, then the long-time dynamics of 
the system can be described in terms of a local, time- 
independent effective Hamiltonian. In this case, the stan- 
dard classification schemes of time-independent topolog- 
ical phases can be applied^— . Therefore, analogous 
phases to those of topological insulators and supercon- 
ductors can be found in dynamically driven systems. Pro- 
posals such as those of Refsl24T-|26l in which topologically 
non-trivial states are induced dynamically in otherwise 
trivial systems, belong to this class. 

This work opens many avenues for future exploration. 
In particular, it will be interesting to search for exper- 
imental manifestations of quasi-energy winding in arti- 
ficial and natural systems. Using quantum walks, it 
is already possible to realize topological phases in one 
dimension^, and we expect that more complex exam- 
ples will soon be accessible. Additionally, the prospect 
of dynamically inducing topological states in naturally 
occuring, topologically trivial, materials is particularly 
intriguing. 

In this paper, we have focused on single-particle dy- 
namics in non-interacting systems. As shown in the ap- 
pendices, the results can also be applied to the adiabatic 
dynamics of gapped many-body systems. It is interesting 
to extend this work to more general many-body systems, 
and to open systems subject to decoherence. These ex- 
tensions could pave the way to finding new methods of 
robust quantum control in many-body systems. 
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Appendix A: Topological classification in the 
presence of disorder 

Topological properties of physical systems are expected 
to be robust against a broad range of perturbations. In 
sections IlII Al and IIII B[ we studied the topological prop- 
erties of evolution operators characterized in terms of ho- 
motopy groups. Non-trivial topology in this homotopy 
group sense is associated with the "winding" of quasi- 



energy in the Brillouin zone. In all the cases considered 
so far, the topological invariants were defined under the 
conditions of (discrete) translational invariance. In par- 
ticular, Eq.® and Eq.([S]) are written in terms of the 
conserved crystal momentum which is associated with 
this symmetry. Below, we generalize the invariant V\ to 
non-translationally invariant situations, using the idea of 
twisted boundary conditions^. The value of the gener- 
alized invariant v\, defined in Eq. (|A2[) below, gives the 
same value as V\, Eq.@, in the absence of disorder. An 
analogous generalization of the invariant 1/3 to the disor- 
dered case can be achieved in a similar manner. 

In the following, we consider a weakly-disordered one 
dimensional system of a finite length La, where a is 
the lattice constant of the "unperturbed" translationally- 
invariant system in the absence of disorder. Disorder is 
weak in the sense that the energy eigenvalues of the in- 
stantaneous Hamiltonian H{t) at each time t remain sep- 
arated into bands which can be associated with the bands 
of the unperturbed Hamiltonian. We assume, as before, 
that the evolution over one full period T only mixes a fi- 
nite number, m, of these bands of the initial Hamiltonian, 
-ff(O), and seek to characterize the topological properties 
of the finite dimensional mL x mL evolution operator 
U(T). 

To probe the topological properties of U(T), we im- 
pose "twisted boundary conditions" parameterized by 
an angle 9, such that all wavefunctions must satisfy 
(x = 0\ip) = e l6 (x = L\ijj). In principle, for each value 
of one can find the corresponding twisted boundary 
condition eigenstates of the Floquet operator U(T), and 
their associated quasi- energies. However, the inconve- 
nience of dealing with twisted boundary conditions can 
be eliminated by studying the transformed Floquet op- 
erator 

U e (T) = e~ lie/L U(T)e l£e/L , (Al) 

where Ug{T) acts on states with periodic boundary con- 
ditions. 

In terms of the twist angle 9, we define the generalized 
topological invariant v\ as 

vi = ^ J d9Tv [Ue(T)-Hd e U e (T)] , (A2) 

where the trace is taken over all states in the mL- 
dimensional low-energy subspace (with periodic bound- 
ary conditions). Notice that expression (|A2j) for v\ is 
analogous to that of v\, Eq.j2|), with the phase twist 9/L 
playing the role of the crystal momentum. In this way, 
the entire disordered system of length L plays the role 
of one giant unit cell of a much larger, periodic system 
with mL bands. As above, v\ counts the net winding of 
quasi-energy as 9 is taken from to 2ir (see Fig[S]). 

Similarly, the analogous generalization of ^3 for a finite 
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tt/T 




-k/T 



FIG. 8: The generalized topological invariant u\, see Eq. 
(IA2|) , counts the net winding of quasi-energy in one dimen- 
sion as the twist-angle 6 is taken from to 2tt. Here, for each 
0, we plot the mL quasi-energy eigenvalues of the evolution 
operator. 



system in three dimensions is given by 

d 3 e 



"3 = 



24tt 2 



e aM x 



(A3) 



Tr 



{U^deJJe) {U^dopUo) (U^dg^Ug) 



where a, /3, 7 label the Cartesian directions {x, y, z}, 6 — 
(9 X ,9y,6 z ) is a vector of twist angles for the boundary 
conditions in each of the three spatial directions, and the 
integrals are taken over — tt < 6 a < tt. 

The generalized topological invariants v\ and £3, de- 
fined in terms of twisted boundary conditions, can also 
be applied to many-body systems, as long as there is a 
finite gap between the ground state and the rest of the 
spectrum throughout the driving period. In particular, 
it can be applied in the case of band insulators. In these 
cases, v\ and £3 are meaningful even in the presence of 
many-body interactions^. 



Appendix B: Relation between quantized charge 
pumping and quasi-energy winding 

In this Appendix, we study the topological struc- 
ture of quantized adiabatic pumping, first proposed by 
Thoulessi^, as described in terms of the winding number 
v\. Consider spinless fermions moving in a periodically 
varying, one-dimensional lattice potential. We assume 
that, at any time t, the m lowest bands of the instanta- 
neous Hamiltonian H (t) are filled, i.e. the system forms 
a band insulator. This is the case if the rate of change 
of the potential is slow enough such that the evolution is 
adiabatic with respect to the band gap at all times, and 
no particles-hole excitations created. Thouless proved 
that, under these circumstances, the number of particles 
transported through the system over one driving period is 
quantized as an integer—. Here, we show that this quan- 
tized transport can be naturally understood as a conse- 



quence of the winding of quasi-energies of the correspond- 
ing evolution operator, as captured by the v\ invariant. 
In fact, such quantized charge transport amounts to a 
realization of the quantized average velocity (or displace- 
ment) discussed in Section UlI Al 

We begin by relating v\ to the charge current inte- 
grated over one period of the dynamics. Assuming that 
the particles do not interact, we study the properties of 
the single-particle dynamics. For simplicity we assume 
translational invariance, but for interacting or disordered 
systems, pumping of the many-body state can be studied 
using the twisted boundary condition method described 
in Appendix [A"l 

It is useful to work in the basis of cell periodic states 
\u k ,a) = e~' lkx \k, a), where \k,a) is the Bloch wave 
function with crystal momentum k and band index a. 
The corresponding transformed Hamiltonian is given by 
Hk(t) — e~~ lkx H (t)e lkx , and the evolution operator is 
given by U k {t) = e~ lkx U{t)e lkx . The Floquet evolu- 
tion operator Uk(T), which acts on the wave functions 
{|ufc,a)} with periodic boundary conditions, only mixes 
the filled bands due to the adiabaticity condition of the 
dynamics. Therefore, when restricted to the manifold of 
filled bands, and projected onto the subspace with a par- 
ticular crystal momentum fc, Uk(T) reduces to U k {T) as 
defined in the main text. 

The instantaneous pumped current in the many-body 
insulating state can be expressed as 



J{t) 



E 
E 

It^OCC 

E 



afEocc. 



dk 
2^ 



[Uk,a\ \Uk,a) 



Ot 



l i[H k {t),x}U k {t)\u k ,a) 



^(u k ,a\U k (ty 
Zir 

dk ~ / 

— (tifc.alC^t)- 1 (d k H k (t))U k (t)\u k ,a) 



where we have used the Heisenberg representation, and 
the summation of a is over the occupied bands. Using the 

relation d k U k (T) = - ^ U k (T - t) (id k H k (t)^j U k (t)dt, 

the total pumped current integrated over one period is 

i-T 

dtJ(t) = (Bl) 



E 

a^occ. 



— (u k ,a\U k {T)- 1 

Z7T 



id k U k (T))\u k ,a). 



Equation (|B1[) is equivalent to the expression of v\ in 
Eq.©, with the understanding that the trace is taken 
only over the occupied bands. As pointed out in Sec. 
IIII Al v\ counts the winding number of the quasi-energy 
bands; hence the total pumped charge can take only in- 
teger values. 

To demonstrate the quasi-energy winding of a system 
that displays quantized pumping, we consider the follow- 
ing example: 



Hit) 



2m 



+ Acos(27r£ — ujt), 



(B2) 
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where we have taken the lattice constant to be 1. The 
cosine potential is used for concreteness, although the 
following argument applies to any moving periodic po- 
tential V(2ttx — Lot). The (real) momentum is denoted 
by p, which should not be confused with crystal momen- 
tum k. 

The evolution operator under the time-dependent 
Hamiltonian above can be exactly obtained by a Galilean 
transformation to the comoving frame. It is straightfor- 
ward to check that the following evolution operator sat- 
isfies the equation of motion id t U(t) = H(t)U(t); 



U(t) 



g — iwtp/ (27r) -iH t 



(B3) 



Here, Hq is the static Hamiltonian Hq — f— — %2 + 
Acos(27rx). The evolution operator after one period can 
be written as 



U(T) 



-iH T 



(B4) 



where T — 2ir/u>. Let {| £;,«)} be the Bloch eigenstates 
for Hq, where k is the crystal momentum and a is the 
band index. By definition, \k,a) is an eigenstate of the 
shift operator e~ 1 ^ with eigenvalue e~ lk . Therefore, we 
may write 



U k {T)\k,a) =e~ lk e 



(B5) 



where Ef., a is the energy of the Bloch state \k, a). Since 
Ek,a is a periodic function of k, it is clear from this ex- 
pression that each quasi-energy band of U(T) is charac- 
terized by a winding number of 1. 

Note that in the example above, the adiabaticity condi- 
tion was not used - the quasi-energy winding is a straight- 
forward consequence of the Galileian invariance of the 
system. Even if Galileian invariance-breaking perturba- 
tions, such as a small time-independent potential, are 
added, as long as the adiabaticity condition is imposed 
(to -C AE g , where AE g is a typical band gap), particles 
cannot be excited from one band to the other during the 
dynamics. Therefore, the quasi-energy winding of each 
band cannot change. 



Appendix C: Relation between the topological 
invariants of evolution operators and Chern numbers 

The Chern numbers have been used intensively to clas- 
sify insulators and superconductors. The first Chern 
number is defined in two dimensions, and was used to 
explain the quantization of the Hall coefficients in inte- 
ger quantum Hall systems^. More recently, Qi et. alr^ 
have shown that topological invariants of two and three 
dimensional topological insulators are closely related to 
the second Chern number in four dimensions. 

In this paper, we have shown that periodically driven 
systems can be characterized by the invariants v\ and 
j/3, which are written entirely in terms of evolution oper- 
ators. Through dimensional reduction^, we establish a 



connection between v\ and v^, and the first and second 
Chern numbers in 2 and 4 dimensions, respectively. 

Qi et. al*£ classified a family of 2n — 1 dimensional 
insulators indexed by a compact parameter through 
Chern numbers defined in 2n dimensions by regarding 
the parameter as an additional component of crystal 
momentum. In this way, they have shown that the first 
and second Chern numbers defined for a family of 1 and 
3 dimensional insulators are given by 



C 2 

c 4 



1 

27 



dPi{6) 



16tt 2 



80 

OP 3 (0) 
80 ' 



where 



Px(0) 



dkTr[a x (8,k)} 



P 3 (0) = /d 3 k£^Tr (fij-^Kaj^-ak 



and 



a? p (0,k) = H)K,0,a|%K,0,/3) 
£f(0,k) = 8 ki af 



^a" 13 + i[a i ,a j ) aP , 



(CI) 
(C2) 

(C3) 



(C4) 
(C5) 



(C6) 
(C7) 



where e % ^ k is the anti-symmetric tensor and the k in- 
tegration extends over the first Brillouin zone. Here 
{|uk, 0, a)} is the set of filled states of the system, which 
satisfy |uk+G,#, a) — |uk, 6, a), where G is a reciprocal 
lattice vector. Note that the Berry connection a, is a ma- 
trix indexed by the internal state labels, a and j3. The 
trace is taken over the filled bands. 

For periodically driven systems, we interpret the para- 
metric dependence of Bloch states on as coming from 
the evolution of the states under the periodic Hamilto- 
nian. Then, it is possible to write the Chern numbers 
in Eq. (|Cl[) and Eq. (|C2[) entirely in terms of the evolu- 
tion operator of Bloch wave functions after one period, 
Uk(T). We identify the cyclic parameter with the time 
variable t, and write the Chern numbers as 



C2 = ^(Pi(T)-Pr(O)) 



C A = 



16tp 



(P 3 (T) - P 3 (0)) 



(C8) 
(C9) 



Therefore, the Chern numbers have the meaning of the 
change of the "polarizations" Pi and P3 after one period. 

The polarizations Pi(T) and Pa(T) at time t = T are 
related to Pi(0) and Ps(0) at time t = through the evo- 
lution of the Bloch states. If the Hamiltonian changes 
slowly compared to the band gap, it does not mix the 
occupied states with states above the gap. Each occu- 
pied state \uk, 0, a) at time t = evolves into a typically 
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different occupied state \v,k,T,a) at time t — T, given 
by \u k ,T,a) = C/ k (T)K,0,a) = ^ U? a (T)\u k , 0, /?), 
where Uk(T) is the evolution operator for crystal mo- 
mentum k and u£ a (T) = (u k , 0, /3\U k (T)\u w , 0, a). The 
evolution of the Bloch states is characterized by a non- 
abelian Berry connection a"^(T), see Eq.fC6j with re- 
placed by T: 

af (T) = af (0) - i{u k , 0, alE^T)- 1 ^^)!^, 0, (3). 

Substituting this relation into the definitions of the two 
polarizations, it is now straightforward to show that the 
Chcrn numbers in Eq. (|C8[) and Eq. (|C9|) are expressed in 
terms of evolution operators Uk(T) and give the topolog- 
ical invariants of evolution operators vi and in Eq. © 
and Eq.©. 

This construction shows the connections between 
driven dynamics in 2n — 1 dimensions and nth Chern 
numbers for n = 1,2. We emphasize that the dynamics 
is not necessarily adiabatic with respect to the separation 
between the low-lying bands. The only condition of the 
applicability of ^1,^3 is that the dynamics is adiabatic 
with respect to a band gap which separates the low-lying 
bands from all higher states (Fig. [TJ), and therefore only 
the low-lying bands are mixed after one period. 



where Q is a unitary operator and to is an arbitrary time, 
then the evolution operator of the system over one period 
is time reversal symmetric [in the sense of Eq.©] and so 
is the effective Hamiltonian, Eq.©. To show this, we 
write the evolution operator in the form 



U(T,0)= lim e -*A*H(r) e -iA«r(T-At) -i>e -iA«r(o) j 

(D2) 

where U(t\,t2) is the evolution operator from time t\ to 
t 2 , and At = T/N. Using Eq. (|Dl|) . and the fact that 
H(t + T) = H(t), one can show that 



QU(T,0)*& = U(2t + T,2to)\ (D3) 



Then, defining Q = U(0,2t )Q, we get that 



Appendix D: Time reversal symmetry in 

periodically driven systems QU(T 0)*Q^ = U(T 0)^ (D4) 

In this Appendix, we prove that if the time-dependent 
Hamiltonian H (t) of a periodically driven system satisfies 

QH(t + t )*& = H(-t + t ), (Dl) which is Eq.©. 
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